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Abstract 

An updated next-to-leading order (NLO) QCD analysis of the spin asymmetries A±(x, Q 2 ) 
and parton distributions Sf(x, Q 2 ) in longitudinally polarized deep-inelastic lepton-nucleon 
scattering is presented within the framework of the radiative parton model taking into 
account recent experimental results. The theoretical framework and the main features 
of the radiative parton model analysis are briefly reviewed. The small- x behaviour of 
the polarized structure function g^(x,Q 2 ) as well as the shape of the polarized gluon 
distribution 5g(x, Q 2 ) are shown to be still hardly constrained by present fixed target data. 



1 Introduction 

The past year has seen much progress in our knowledge about the nucleons' spin structure due 
to new experimental results on the spin asymmetry A±(x, Q 2 ) ~ g^(x, Q 2 )/F^(x, Q 2 ) in deep- 
inelastic scattering with longitudinally polarized lepton beams off nucleon targets (N = p, n, d). 
In particular, previous sparse and not very precise experimental information on the neutron 
asymmetry A\ |l] has been succeeded by more accurate data by the El 54 collaboration [0 
where also the kinematical coverage in x and Q 2 was slightly extended. Recently the HERMES 
group also has presented first results on A™ || and preliminary data for a proton target (Af) 
have been reported on this workshop [|J. SMC has released a new detailed analysis of A\ |j] 
which does not indicate a rise of g\ at small- a; anymore and, finally, first very accurate, but still 
preliminary results from E155 on g\ were also presented on this workshop 0. 

In view of these recent experimental developments it seems to be worthwhile to reanalyse 
these dataQ in terms of polarized parton distributions 5f in the framework of perturbative 
QCD. Preceding studies H, [K| have revealed that the detailed x-shape of the polarized gluon 
distribution 5g(x,Q 2 ) was only weakly constrained by the data of that time even though a 
tendency towards a sizeable positive total gluon polarization, J 5g(x, Q 2 = 4 GeV 2 )<ix > 1, was 



found H, P, 10]. It it thus interesting to study to what extent these old results are confirmed by 



the new data sets and whether one can further pin down the polarized gluon distribution. 



t Present address 

: The preliminary E155 Q and HERMES Q proton data are not available yet and hence not included in our 
analysis so far. 
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In the following we will concentrate exclusively on calculations of A± to NLO accuracy 
which became possible only recently after the derivation of the required spin-dependent two- 
loop anomalous dimensions |~1], [12| . A first such complete and consistent NLO study has been 
presented in || (based on a corresponding LO analysis [0), where the underlying concept has 
been the radiative generation of parton distributions from a valence-like structure at some low 



bound-state like resolution scale /i. In the unpolarized case this had previously led [14], [15i , e.g 



to the successful prediction of the small-x rise of the proton structure function as later on 
observed at HERA. Other NLO analyses of polarized DIS data can be found in || [l(| . 

In the next section the basic theoretical framework for polarized DIS beyond the leading 
order is briefly discussed hereby defining our notations, and the main features of the radiative 
parton model analysis || are reviewed. Section 3 contains the presentation and discussion of the 
new quantitative NLO QCD results and, finally, the main findings are summarized in section 4. 



2 Theoretical Framework 

Measurements of polarized deep inelastic lepton nucleon scattering yield direct information [1- 
6,16] only on the spin-asymmetry 

a n (t n 2 \ ~ dii^Q^) 9i(x,Q 2 ) m 

1 VM>- Ff(x, Q*) Ff(x, Q 2 )/ [2x(l + R N (x, Q 2 ))} ' 1 } 

where N = p,n,d (in the latter case we have used gf = \{gi + PiXl — |^d) with u) D = 0.058) 
and R = Fl/2xF\ = (F 2 — 2xF\) j1xF\. In NLO, the polarized structure function gi(x,Q 2 ) in 
([[]) is related to the spin-dependent quark and gluon distributions (Sf N ) in the following way: 

g?{*> Q 2 ) = \H Q 2 ) + Q 2 ) + ® + *f) + 5C 9 ® S 9}\ (2) 

q ^ > 

with the convolutions (®) being defined as usual. The NLO pieces entering i.e., 5f N , SC q , 
5C g , depend on the factorization convention (scheme) adopted. Since the two-loop anomalous 
dimensions of fill , H3 refer to the conventional MS dimensional regularization prescription we 
prefer to work also in this scheme. The appropriate spin-dependent MS Wilson coefficients 5C q 
and SC g can be found, e.g., in |IT[ |J. It should be recalled that the first moment T^(Q 2 ) = 
Jo 9i{ x i Q 2 )dx of eq. (0) is in the MS scheme simply given by 



(3) 



where we have introduced the first moments Af N (Q 2 ) of the polarized distributions 6f (x, Q 2 ) 
by defining as usual 

Af N (Q 2 )= f Sf N (x,Q 2 )dx (4) 



where / = u,u,d,d,s,s, and g and have used (see, e.g., [0, ||) 5C q (x)dx = —3Cf/2 and 
Jq 1 5C g (x)dx = 0. Thus, the total gluon helicity Ag(Q 2 ) does not directly couple to T^(Q 2 ) due to 
the vanishing of the integrated gluonic coefficient function in the MS factorization scheme. Other 
factorization schemes [[HJ are of course also allowed but afford a proper scheme transformation 
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(which cannot be uniquely fixed) such that the physical quantity g^ remains scheme independent 
up to 0(a s ). 

The NLO Q 2 -evolution of the polarized parton distributions 5f(x,Q 2 ) (henceforth we shall, 
as always, use the notation 5q p = 5q and q p = q) governed by the two-loop anomalous dimen- 
sions |Tl| , [12| is performed most conveniently in the Mellin-n moment space where the solutions 
of the evolution equations (see, e.g., refs. Hl4|, ||) can be obtained analytically, once the boundary 
conditions at some Q 2 = Qq, i.e., input densities 8f(x,Q 2 ) ) to be discussed below, are specified. 
Furthermore, in Mellin-n space the convolutions in (|2|) reduce to simple products. Having ob- 
tained the analytic NLO solutions for the moments of parton densities or the analogous n-space 
expression of (§) the desired x-space results for Sf(x, Q 2 ) or g^{x, Q 2 ) are then simply obtained 
by a standard numerical Mellin inversion as described, e.g., in fT4}| . 

In our analysis only quarks with m q < Aq C d, i-e., u, d, s, will be treated as light partons in 
the evolution equations while the charm contribution to (as well as to F±) is calculated via 
the appropriate massive 7* g — > cc fusion process [17] although it turns out (see later) that in the 



kinematical region covered by present fixed target experiments [1-6,16] the charm contribution 
is extremely small and practically irrelevant. 

To finally fix the NLO input parton densities Sf(x, Qq) we perform fits only to the directly 
measured spin asymmetry A±(x, Q 2 ) in (|I|), rather than to the derived g^(x,Q 2 ). The main 
reason for that is that in some older experimental analyses gi(x, Q 2 ) has been extracted under 
the assumption of the Q 2 -independence of Af(x, Q 2 ), which is - although presently available 
data do not exhibit any significant Q 2 -dependence within the experimental errors - theoretically 
not warranted due to the different Q 2 -evolutions of the numerator and denominator in (|T|). 

As already mentioned in the introduction, the other main ingredient of our NLO analysis M 



is that we follow the radiative (dynamical) concept [14], y| by choosing the same low input scale 
Qq = fi 2 = 0.34 GeV 2 and implementing the fundamental positivity requirement 

\5f(x,Q 2 )\<f(x,Q 2 ) (5) 

down to Q 2 = [j 2 . The analysis affords some well established set of unpolarized NLO parton 
distributions f(x, Q 2 ) for calculating F^(x, Q 2 ) in (|l|) and as reference distributions in which 



will be adopted from ref. |T5|. 



In addition to (H), the first moments Af(Q 2 ) of the NLO polarized parton distributions are 
taken to be subject to two very different sets of theoretical constraints (IB], |8| related to two 
different views concerning the flavor SU(3)f symmetry properties of hyperon /3-decays. One 
set ('standard' scenario) assumes an unbroken SU(3)f symmetry between the relevant matrix 
elements leading to the following sum rule constraints, 

Ag 3 = Au + Au- Ad- Ad = g A = F + D = 1.2573 ±0.0028 (6) 
Aq 8 = Au + Au + Ad + Ad- 2(As + As) = 3F — D = 0.579 ± 0.025 (7) 



with the values of g A and 3F — D taken from |TB[ . It should be noted that the flavor non-singlet 



combinations Aqs t s in (|6|) and (0) remain Q 2 -independent also in NLO JIT, ^2, 

As a plausible alternative to the full SU(3) / symmetry between charged weak and neutral 
axial currents required for deriving the 'standard' constraints ([]) and (|^), we consider a 'valence' 



scenario [13, |T!| where this flavor symmetry is broken and which is based on the assumption [19 



that the flavor changing hyperon /5-decay data fix only the total helicity of valence quarks: 

Au v (fi 2 ) - Ad v (fi 2 ) = g A = F + D = 1.2573 ± 0.0028 (8) 
Au v (fi 2 ) + Arf y (/i 2 ) = 3F — D = 0.579 ± 0.025 . (9) 
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Figure 1: Comparison of our reanalysed NLO 'standard' results for Af(x,Q 2 ) (solid lines) with all 
presently available data [1-3,5,16]. The Q 2 values adopted here correspond to the different values 
quoted in [1-3,5,16]. Also shown are the results of our previous analysis @ (dashed lines). 



The 'standard' scenario always requires a finite total strange sea helicity of As = As ~ —0.05 
in order to account for the experimentally observed reduction of with respect to the Gourdin 
and Ellis and Jaffe estimate |2(J . Within the 'valence' scenario, on the contrary, a negative light 
sea helicity Au = Ad = Aq ~ —0.07 alone suffices and we shall assume a maximally SU(3)j 
broken polarized strange sea input 5s(x,fi 2 ) = 5s(x,fi 2 ) = here M. Finally, we note that in 
both above scenarios the Bjorken sum rule manifestly holds due to the constraints (|6]), (|8|). 

3 Quantitative NLO Results 

Turning to the determination of the polarized NLO parton distributions Sf(x, Q 2 ) it is helpful to 
consider some reasonable and simple, but still sufficiently flexible ansatz for the Sf(x,fj, 2 ) with 
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Figure 2: The polarized NLO MS densities at Q 2 = 4GeV 2 in the 'standard' scenario as obtained in 
our new and old || analyses. Also shown are the distributions obtained in two other fits employing 
some extra constraints on the input distributions. It should be noted that as a result of the fits the 
sea always turns out to be SU(3)f symmetric, i.e., q = u = d = s. 

not too many free parameters. As a general ansatz we take || 

5/(x,/i 2 ) = N f x a '(l-x)P'f(x,ti 2 ) (10) 

where / = Uy, dy, q = u = d, s, g which links the polarized input distributions to the unpolarized 
ones as taken from ]15] (in this way the positivity requirements (|J) can be trivially implemented). 
One can argue that the available inclusive DIS data do not allow for a fully flavor-decomposed 
ansatz like (|T0|) but instead only for a separation into non-singlet, singlet and gluon distributions. 
This is of course true, but the aim of our analysis is not only to simply fit the available data but 
also to provide a realistic set of parton distributions which can be further probed in processes 
other than DIS, i.e., to have 'predictive power'. For obvious reasons, we have not taken into 
account any SU(2)f breaking input (5u ^ 5d) in (JITj). Moreover the number of parameters in 
( |10"D can be further reduced without any worsening of the fit by setting (3 Uv = j3d v = 0, a s = 0, 
and (3 S = in flip]). 

A comparison of our best fit employing the 'standard' scenario constraints @ and (0) with 
the available data on A±(x,Q 2 ) [1-3,5,16] is presented in fig. 1. Also shown are the results 
for obtained by using our previous fit results [|J. The results in the 'valence' scenario are 
indistinguishable from the ones shown and hence suppressed. The total x 2 °f the 'new' and 
'old' optimal 'standard' scenario fits is 123.02 and 144.37, respectively, for 168 data points. As 
can be inferred from fig. 1 the results for the proton (A^) as well as for the deuteron (Af) are 
basically not affected by the reanalysis despite of the rather large differences in the neutron case 
(A™) mainly due to the rather precise, new E154 data 0. This can be understood better by 
comparing the individual parton distributions 5f(x,Q 2 ) rather than A^ itself. This is done in 
fig. 2 for the 'new' and 'old' NLO MS densities at Q 2 = 4GeV 2 . As can be seen from fig. 2 
apart from a slightly smaller gluon distribution Sg only Sdy changes by quite a large amount. 
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Figure 3: The small-a: behaviour of g\ (upper part) and <?" (lower part) in NLO at Q 2 = 10 GeV 2 as 
predicted by various 'standard' scenario fits employing different boundary conditions for the polarized 
input gluon distribution Sg(x, fx 2 ). For the optimal fit input also the effect of not including the charm 
contribution to g^ n is shown. 

This is because 5dy is mainly probed in g™ (i.e., A™) where we have observed the largest changes 
as compared to our old results || in fig. 1. 

Also shown in fig. 2 are the results of two other fits which are based on some additional 
constraints on the input distributions. For the '(3f = 0' fit we have set (3f = in our ansatz 
(|T0"D. The total x 2 °f 123.6 is very similar to our best fit and this confirms the results of a 



recent fit by the E154 collaboration |21] based on our radiative parton model framework |y 
where all /3's where chosen to be zero from the very beginning. The 'static Ag' fit also yields 
a similar total \ 2 (124.28) and the idea behind this fit deserves a further explanation. If one 
studies the Q 2 -evolution of the total gluon polarization Ag(Q 2 ) one observes that Ag rises for 
increasing values of Q 2 if one starts with an input Ag(fi 2 ) > Ag static . On the other hand Ag 
decreases for increasing Q 2 if Ag(fj, 2 ) < Ag static hence remaining constant for all values of Q 2 for 
Ag{ji 2 ) = Ag static (cf. fig. 5). In LO the precise value of Ag static can be easily obtained from the 
evolution equation for Ag by setting dAg(Q 2 )/d\nQ 2 = 0. This yields 



A<7. 



static 



— AS 
2 



O(-0.15) 



where AS = ^ (Aq + Aq) is the total helicity carried by quarks and antiquarks and (TlT|) is only 
subject to a small correction in NLO. 
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Figure 4: The experimentally allowed range of NLO polarized gluon densities at Q 2 = 10 GeV 2 for 



the 'standard' scenario. Also shown is the unpolarized gluon distribution of ref. [15 



Fig. 2 also reveals that the valence (and to some extent also the sea) distributions are quite 
well determined by present data in our new fits despite of the different underlying constraints for 
the input distributions. On the contrary the polarized gluon density 5g is still hardly constrained 
at all, even the rather small and 'exotic' Sg resulting from the 'static Ag' fit is not excluded. 

Inevitably the large uncertainty in 5g implies that also the small-x behaviour of g\ beyond the 
experimentally accessible x-range is completely uncertain and not predictable as is demonstrated 
in fig. 3 for the proton and the neutron case in the 'standard' scenario for Q 2 = 10 GeV 2 . Apart 
from the optimal input and the already discussed 'static Ag J fit we also present in fig. 3 the 
results obtained by chosing three other extreme boundary conditions for Sg(x,n 2 ). The total 
X 2 values for these '5g = -g% '5g = 0', and '5g = g 7 inputs are 134.68, 124.24, and 127.44, 
respectively. Even these inputs give still excellent fits and only the largest possible negative 
input in the radiative parton model (cf. eq.(|5|)), i.e., 5g(x,{i 2 ) = —g(x,fi 2 ), is disfavored by its 
X 2 value. As can be seen, all these extreme inputs give indistinguishable results for gf' n in the 
experimentally covered x-region (x > 0.01) whereas they lead to a rather large spread in the 
small- x region making any predictions impossible here. Note that one obtains completely similar 
results also in the 'valence scenario'. This uncertainty implies also a large theoretical error from 
the extrapolation x — > when calculating the first moments T^ ,n . A conservative theoretical 
estimate for rf' n , taking into account the maximally allowed spread in g\' n by the positivity 
requirement (B) in the radiative parton model, i.e., '5g = —g' . . .' 8g — g' inputs, yields 



T\(Q 2 = 10 GeV 2 ) = 0.133 ± 0.008 , r?(Q 2 = 10 GeV 2 



-0.062 ±0.008 



(12) 



Also shown in fig. 3 is the effect of not including the charm contribution to gf n for our 
optimal fit results. As already mentioned in section 2, g^ harm as calculated from the appropriate 
massive 7*g — > cc subprocess |17| is negligibly small in the experimentally covered x-region and 
remains small even for lower values of x. 
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Figure 5: NLO (Revolution of the first moments Ag, AS, and r^' n as well as of the total orbital 
angular momentum contribution L z to the helicity sum rule fll3p from our low input scale /i 2 up to 
Q 2 = 1000 GeV 2 for our 'old' [Q and 'new' optimal fits as well as for the 'static Ag' input. Also 
shown are some recent experimental results for and from SMC |5[] and E154 0, respectively. 



Fig. 4 compares the different NLO gluon distributions at Q 2 = 10 GeV 2 as obtained from the 
various inputs discussed above. It is immediately obvious that other measurements apart from 
DIS are required to further pin down 5g. For completeness the total polarizations Ag(Q 2 ) at 
Q 2 = 10 GeV 2 for the '6g = g\ best fit 6g, '6g = 0', and 'static Ag' inputs are 3.2, 1.45, 0.31, 
and -0.12, respectively. The value for the best fit gluon input is quite close to the value of 1.74 
found in our previous analysis [|| and is in agreement with other fits || [10], pl|p|. 

Finally, let us turn to the NLO (Revolution of the first moments Ag(Q 2 ), AS(Q 2 ), and 
r^' n (Q 2 ) which is shown for the 'standard' scenario in fig. 5 for Q 2 values ranging from our low 
input scale ft 2 up to Q 2 = 1000 GeV 2 . Also shown is the Q 2 -dependence of the total orbital 

2 It should be noted that the possibility of having a slightly negative Ag was recently also found in [2^] but 
without discussing the special 'static' properties of such a boundary condition. 
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angular momentum L z as obtained from the helicity sum rule 

l -= l -/^{Q 2 ) + Ag{Q 2 ) + L z {Q 2 ) . (13) 

It is interesting to observe that at our low input scale for our optimal fit (this hold true also for 
our previous results M) the nucleon's spin is dominantly carried just by the total helicities of 
quarks and gluons, i.e., L z (fi 2 ) rj 0, and only during the (Revolution a large negative L Z (Q 2 ) is 
being built up in order to compensate for the strong rise of Ag(Q 2 ) in (|13|) , see fig. 5. 

For our 'static Ag' fit the situation is completely different. First of all one should note that 
Ag(Q 2 ) is indeed independent of Q 2 and quite small (cf. eq. flTTD ) and this in turn implies that 
also L z is practically Q 2 -independent because AE(Q 2 ) is only weakly Q 2 -dependent in the MS 
scheme (it should be recalled that in LO AX(Q 2 ) = const.). But more striking is the fact that for 
this boundary condition the quark and gluon contributions to the helicity sum rule ( |T3"D cancel 
each other (see eq. (|TT|)) implying that for all values of Q 2 the proton spin is entirely of angular 
momentum origin. This result is quite puzzling and completely different from the intuitively 
expected vanishing of L z at some low bound-state-like scale as observed for our optimal fit but 
cannot be excluded yet by the presently available fixed target data. 

4 Summary and Outlook 

We have presented an updated NLO QCD analysis of the DIS spin asymmetry A^(x,Q 2 ) data 
in the MS scheme in the framework of the radiative parton model. Compared to our previous 
results || we have observed a rather large change in the polarized Sdy distribution which is 
mainly due to new data for A", in particular from E154 ||. In contrast to the polarized valence 
(and sea) quarks the gluon density 5g(x, Q 2 ) turns out to be still hardly constrained at all by 
present data. Our optimal fits, however, still favor a rather sizeable total gluon helicity, e.g., 
Ag(Q 2 = lOGeV 2 ) ~ 1.45, but it was shown that even rather exotic boundary conditions for 
Sg, with a rather small total helicity, such as the 'static Ag' input, yield excellent descriptions 
of all available A^(x, Q 2 ) data. 

The latter gluon input has the striking consequence that the spin of the nucleon is entirely 
made of orbital angular momentum for all values of Q 2 contrary to the intuitively expected 
vanishing of L z at some low scale fi 2 as observed for our optimal fits. Finally, it was shown that 
the uncertainty in 5g induces a rather large spread in the small- x behaviour of g^ making any 
predictions impossible. An estimate for the theoretical error in the determination of T^ ,n due to 
the x — > extrapolation uncertainty was given. 

Future fixed target DIS data, in particular from E155, as well as semi-inclusive measurements 
from SMC and HERMES will help to pin down the polarized quark distributions more precisely 
but it cannot be expected that 8g can be further constrained by such measurements since the 
lever-arm in x and Q 2 is too limited for an indirect determination of Sg from scaling-violations. 
A realization of the currently discussed upgrade of HERA to a polarized ep collider would be 
extremely helpful in this respect to gain more insight into Sg as well as to pin down the small-x 
behaviour of g{(x, Q 2 ) more precisely. 
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